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Abstract 

Drawing inspiration from flight behavior in biological settings (e.g. territorial bat¬ 
tles in dragonflies, and flocking in starlings), this paper demonstrates two strategies 
for coverage and flocking. Using earlier theoretical studies on mutual motion camou¬ 
flage, an appropriate steering control law for area coverage has been implemented in a 
laboratory test-bed equipped with wheeled mobile robots and a Vicon high speed mo¬ 
tion capture system. The same test-bed is also used to demonstrate another strategy 
(based on local information), termed topological velocity alignment, which serves to 
make agents move in the same direction. The present work illustrates the applicability 
of biological inspiration in the design of multi-agent robotic collectives. 

Index Terms - Path Planning for Multiple Mobile Robots or Agents, Autonomous Naviga¬ 
tion, Distributed Robot Systems, Wheeled Robots 

1 Introduction 

The last few decades have witnessed an increase in research efforts towards uncovering mech¬ 
anisms behind pursuit behavior and collective motion IM in nature. Parallel to this 
stream of endeavors, a number of mathematical models and appropriate feedback mecha¬ 
nisms [TnUTS] have been introduced to bring these ideas from natural settings to unmanned 
vehicular systems. Moreover, some groups in the robotics community have performed suc¬ 
cessful implementation of these concepts EnHEQ], and thereby demonstrated the power of a 
bio-inspired approach towards synthesizing collective motion. This current work of ours is 
similar in spirit, and provides an indoor demonstration of two strategies for coverage and 
flocking behavior. 

*This paper is an extended version of ,Tj. 
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Figure 1: Mobile robot based experimental platform (Pioneer 3 DX) with two-wheel differ¬ 
ential and caster. 


In this paper, we have implemented mutual motion camouflage (MMC) [13] on a mobile 
robot test-bed. Existing literature on dragonflies [2T] provides qualitative analysis of territo¬ 
rial battles, wherein the trajectories display spiraling motion consistent with the theoretical 
predictions [22]. This particular bio-inspired control algorithm inherits an appealing cover¬ 
age property through the mechanism of space filling curves, and our implementations are 
able to reproduce coverage patterns similar to the predicted ones. 

Although there has been a long history of control algorithms for flocking, almost every 
model of collective motion predicts diffusive transport of information. But, contrary to the 
existing models, recent findings [18] from starling flocks suggest that directional information 
within a flock propagates with an almost constant speed, and this linear growth of infor¬ 
mation can be explained by models with wave-like aspects. In the later part of this paper, 
we have introduced a control strategy, called topological velocity alignment (TVA), which 
conforms to this criterion and can explain how information about local neighbors can influ¬ 
ence the agents in a flock to align their headings in a single common direction. Moreover, 
an ongoing research work has shown that empirically observed curvature values in starling 
flocks carry a strong resemblance with the curvature values predicted by the steering control 
laws associated with TVA. Hence it seems reasonable to use TVA strategy for collective 
motion synthesis. Furthermore, our implementation results in real robots have shown that 
reduction in neighborhood size and external perturbation (similar to predator attack) can 
split a flock into smaller subgroups. 

2 Self-steering Particle Model 

By considering an agent to be a unit-mass, self-steering particle, we use natural Frenet frame 
equations [23] to describe its motion. This approach to describe an agent’s motion is based 
on the fact that, while the unit tangent vector x,; for a given trajectory rj is unique, we can 
choose two unit vectors (y*, z*) such that {x*, y*, z*} defines a right-handed orthogonal frame. 
The evolution of this frame along the agent’s trajectory is given by 

r i = I'jX; 

= Vi + v i*i) 
y i = — r'jUjXj 

Z i = — I'jUjXj, 
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where zy represents speed and ( Ui,Vi ) are the natural curvatures. In a similar way, by re¬ 
stricting the motion of an agent to a planar setting, we can model its motion as 

r, = zyx,; 

Xj = (2) 

y t = -ViUiXi, 

where Xj denotes normalized velocity of the z-th agent and y* is the orthogonal rotation of 
Xj in the counter-clockwise direction. It should be noted here that this way of modeling a 
trajectory requires only twice-differentiability of the trajectory. 

Remark 2.1. As we use Hilare type mobile robots (Fig [Ip as a test-bed for control 
algorithms, ([2]) provides a natural choice to model their motion. 

3 Mutual Motion Camouflage (MMC) 

Here we consider the mutual interaction between two agents each applying the same pursuit 
law, while perceiving the other one as a target. As the dynamics of MMC in a planar 
setting has been studied earlier [13] , we just reiterate some key results in order to have 
a comprehensive framework. Allowing different speeds for the agents, we begin with the 
following symmetry: 

UlUi = U2V-2 = U. (3) 

Then the dynamics of the relative motion vectors, namely r = ri r 2 , g = r = zqxi — Z/ 2 X 2 
and h = g 1 = r , can be expressed as 


r = g 

g = uh (4) 

h = —u g. 

Now we introduce three scalar shape variables defined as p = |r|, 7 = (r • g)/|r| and A = 
(r • h)/|r|. Then, according to TT)]. we have 



where, p > 0 denotes the feedback gain. As shown earlier, the dynamics of relative motion 
flU can be reduced to yield a second order dynamics given by 

p = 7 

7 = ( 1 /p - p) (6 2 - 7 2 ) , 

where, 6 = |g| = |h| is conserved along any trajectory of (J4j). As detailed in the original work 
[IB] , individual trajectories can be reconstructed from the solutions of dHJ). Moreover, the 
solutions of the reduced dynamics ([ 6 ]) constitute level sets for another conserved quantity, 
defined as 

E{P, 7) = p 2 (h 2 - 7 2 )e" 2w = E(p 0 , 70 ). (7) 
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3.1 Stabilization through Dissipation 

As it will be described in section [5l implementation of the original MMC law (J5]) in a labo¬ 
ratory environment did not result in a satisfactory performance. The trajectories diverged 
from the theoretical predictions, and the errors kept building up (Fig [4]) . This type of be¬ 
havior can be attributed to the lack of damping in the reduced dynamics (J6]) . As a result, 
every small error was kept unchecked, and accumulation of these errors gave rise to a poor 
performance of the MMC feedback law. Based on this observation, we considered a modified 
version of the feedback law which introduced a dissipation term. As E(p, 7 ) was not staying 
constant during an implementation of MMC feedback law ([5]), we added a dissipative term 
in the feedback control to counter-act any deviation from the predicted trajectories. The 
resultant control law can be expressed as 


utot = u + u dis = -pX + k d \i{E{p, 7 ) - E d ) , ( 8 ) 

where E d is set as the initial value of the conserved quantity E(p, 7 ) (E d = E(p 0 , 70 )). Earlier 
research | 12 | has shown that the modified control law (jBj) with k d > 0 makes the periodic 
orbit (with energy E d ) orbitally asymptotically stable, and the corresponding domain of 
attraction is characterized by {(p, 7 ) : p > 0 , —8 < 7 < 8, (p, 7 ) A (l//q 0 )}. 

4 Topological Velocity Alignment (TVA) 

Here we formalize the strategy of topological velocity alignment (TVA), and assume that each 
member in a group of n-agents uses this strategy to move together while keeping its heading 
parallel to the neighborhood center of mass velocity. Letting A/j denote the neighborhood of 
the i-th agent, the center of mass (COM) velocity of this neighborhood is given by 

V COM = 7^^ ? A x i> (9) 

1 d jeMi 

where |A/)| represents the number of neighbors influencing the z-th agent. Next, by assuming 
that v COM does not vanish to zero, we define the direction of the center of mass motion as 

X M = |-| ■ (10) 

I V COM I 

It should be noted that x_/y. is not well-defined over a thin set in the state space. As the 
chance of getting into the thin set is essentially zero , we can overlook this situation for all 
practical purposes. Now we introduce a contrast function 

@i = ^ ( X M - Xi) • (x M - Xi) = 1 - Xi ■ x M , ( 11 ) 

as a quantitative measure for the misalignment between the heading of an agent and the 
direction of motion of its neighborhood center of mass. Clearly, this contrast function (©*) 
assumes its minimum value (= 0 ) whenever the z-th agent’s velocity is aligned with its 
neighborhood center of mass velocity, and increases monotonically with increase in the mis¬ 
alignment between them. Thus, 0,; can be interpreted as a measure of departure from our 
goal of achieving alignment. 
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Next, by assuming a non-zero velocity for the neighborhood center of mass (v OOM A 0), 
we propose a control law 


Ui = n 

Vi = pL 


xy, • y i 

Vi 

x A/i ' z * 
Vi 


( 12 ) 


where /i > 0 denotes a positive gain, and y*, z* carry their usual meaning. Alternatively, 
lateral acceleration for this choice of control laws (1T2|) can be expressed as 


a“ = fx,/. - [x_v, • Xj]xj), 


(13) 


and this provides a physical intuition behind (TT2]) as the lateral acceleration is proportional 
to the projection of the normalized velocity of its neighborhood center of mass onto the 
transverse of its own direction of motion. 


Remark 4.1. Earlier works IF7[\25] have considered a very similar form of this control law 
with three components for attraction (while the agents are far away), repulsion (to avoid 
collision) and velocity alignment. However, the control law introduced in this current study 
considers only velocity alignment, but extends the scope from a planar setting to a three 
dimensional environment. Moreover it relaxes the assumption on uniform speed of the col¬ 
lective by allowing the agents non-uniform and time-varying speed profiles. This relaxation 
plays an important role in the context of applying this control law to a group of heterogeneous 
agents. 


4.1 Special Case: A planar 2-agent system 

As analysis of an ?vagent system with neighborhood defined as the set of A'-nearest neighbors 
poses hard challenges, we begin by considering a special case, namely a two-agent system 
with the motion restricted to a planar setting. Clearly, the neighborhood center of mass 
never loses speed in this case (as the neighborhood is nothing but the other agent). 


4.1.1 State space and its reduction onto the shape space 

By modeling the dynamics of an agent using natural Frenet frame equations restricted to a 
planar setting (J2]) , its position and the frame vectors can be packed inside a 3 x 3 matrix 
gi e SE( 2), the Lie group of rigid motions in a plane. By excluding collocation of the agents, 
we define the state space of the system as 


AA state 


g u g 2 eSE(2)xSE{2) 


01 e;s ¥= 02 e 3 


(14) 


where e 3 = [0 0 1] T and g, = 


x* 

0 


y; 

0 


r* 

1 


In terms of this Lie-group formulation the system dynamics can be represented as 


9i = gi&(ui) = g i u i (A 1 u i + A 2 ) (15) 

where A\ and A 2 represent infinitesimal generators of planar rotation and translation, re¬ 
spectively. Moreover A± and A 2 can generate the Lie-algebra under bracketing. As we are 
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interested in steering laws which leave our system dynamics invariant under rigid motion, we 
can formulate a reduction to the shape space , a quotient manifold M. state /SE(2 ) of relative 
positions and velocities of the agents. We define g e SE(2) as 


9 = 9i l 92 


X! • x 2 X! • y 2 X! • r 
yi • x 2 yi ■ y 2 yi ■ r 
0 0 1 


(16) 


where r = r 2 iq denotes the baseline vector. Now the shape space for the planar two-agent 
system can be defined as 


M 


shape 


9 E SE( 2) gf 3 + p 23 ^ 0}. 


(17) 


Moreover, g assumes a left-invariant dynamics on SE(2) as 


9 = 9i (18) 

where £ = £ 2 (w 2 ) — <7 _1 £i(iti)<7 e se(2), and the proposed control law (fT2l) depends only on 
the shape variable g as 


Therefore from B31 it can be concluded that the reduced dynamics (TT8|) evolves on the shape 

Space shape- 



Figure 2: Illustration of scalar shape variables (p, used to parametrize the shape space 

Mshape- 


*6 {1,2}, i = {l, 2 }\{i}. 


(19) 


4.1.2 Shape dynamics 

Now, through polar parametrization, we introduce some geometrically meaningful scalar 
variables to parametrize the shape space. By defining 

r = r 2 — ri = pe 

Xi = e i01 (20) 


we introduce i/j(= n — •& + 6 \) and </>(= 9\ — 62 ) to represent the relative orientation of the 
velocity vectors. Clearly, i/j represents the relative orientation of x t with respect to the 
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baseline vector r, and (p represents the misalignment in velocity directions. From (jT6|) one 
can notice that g e Ai shape can be represented in terms of these scalar shape variables as 


9 = 


COS 0 

— sin 0 


0 


sin 0 

cos 0 


0 


—p cos ip 
p sin ip 
1 


( 21 ) 


Then, by a straightforward calculation, one can show that the shape dynamics are given in 
terms of the scalar variables by 


p = v i cos ip — v 2 cos(0 — <p) 

ip = V\ U\ - \v\ sin ip — v 2 sin pip — 0)1 (22) 

P 

<P = V\U\ — V 2 U 2 . 


4.1.3 Analysis of topological velocity alignment 

Here we consider a particular context of the two-agent planar system wherein each agent 
employs the strategy for topological velocity alignment (TVA). In terms of the original state 
variables the contrast functions take the form 

0i = (x; - Xj) • (xj - Xi) , !6{ 1,2} (23) 

and the z-th agent is declared to attain TVA if (-), = 0. From (l23l) one can immediately 
notice equality of the contrast functions for both agents, and hence we define a common 
contrast function © = ©i = © 2 . Clearly, this common contrast function can be represented 
in terms of scalar shape variables as 


© = 1 — cos (p, 


(24) 


and hence we can conclude that 9 = 0*/ and only if <p = 0. 

Next, for this two-agent planar system we define a 2-dimensional TVA manifold M.tva c 

Mshape 


■M-TVA — | Pi 0, <P e M-shape <P — 0 j 


(25) 


Moreover, from (TTUT) one can notice that the feedback law can be expressed in terms of shape 
variables, taking the form 


Mi = — (FT) sin 


u 2 = 



(26) 


Clearly, the steering control becomes identically zero on the TVA manifold, and as a conse¬ 
quence the mismatch in velocity direction remains identically zero (122]) . Now we will formally 
introduce the notion of invariance. 


Definition 4.2 (Invariant Manifold). A manifold M. is said to be invariant under the flow 
of a vector field X on Xi if for any m e M., F t (m) e M for small t > 0, where F t (-) denotes 
the flow of X. One can show that this condition is equivalent to X being tangent to the 
manifold Ai. 


7 





If both agents employ a steering control of the form ([26]) . the closed loop dynamics for a 
two-agent planar system can be represented as 

fj = U\ cos ip — v 2 cos(ip — 0) 

0 = —p sin 0- \u\ sin 0 — z/ 2 sin(^> — 0)1 (27) 

P 

0 = — 2/rsin0. 


We should note that prohibition on collocation, i.e. p > 0, is not enforced by these dynamics. 
Now we state the following proposition associated with the closed loop shape behavior. 

Proposition 4.3. The topological velocity alignment manifold A4tva c M-shape is invariant 
under the closed loop shape dynamics (1271) . Moreover if'yft) e M. shape is a trajectory of (1271) 
which does not have a finite escape time, and 0(0) A 2, then 

0(f) —> 0 as t -* oo, (28) 


i.e. 7(f) converges asymptotically to M.tva- 
Proof. From (124|) and (127|) we have 

0 = 0 sin 0 = —2p sin 2 0 = —2/i@ (2 — 0) (29) 


and hence M.tva c M-shape is an invariant manifold under the closed loop shape dynamics. 
This equation (1291) also implies that the agents will keep on moving in the opposite direction 
(0(f) = 2) if their initial directions were opposite to each other (0(0) = 2). 

Moreover, by assuming 0(0) e (0,2), we have 


0(f) 


2e" 4 ^ 

C + e~^t 


where the constant C is defined as C = Q^jy — 1. Since e —» 0 as f —> 00, we have 0(f) —> 0 


as f —* 00. 


□ 


Next we focus on the restricted dynamics on Mtva, and by substituting 0 = 0, (127[) 
yields 

p = (o\ ~ v 2 ) cos ip 


•ip = — (v\ — v 2 ) sin?/;. 
P 


(30) 


Now, by assuming zq — z / 2 to be non-zero, the evolution of a phase plane trajectory can be 
described as 


dp p cos ip 
dip sin ip ’ 


(31) 


which can then be integrated to show that the quantity p[t) sin ip[t) is conserved along any 
trajectory on the TVA manifold ( M.tva )• Hence it is clear that the region {p > 0,7r > ip > 0} 
(or similarly {p > 0, 0 > ip > — 1 r}) is positively invariant under the restricted dynamics (130|) . 
The corresponding phase-portrait is shown in Fig [3] 





4.2 Algorithm for an n-agent system in a three-dimensional set¬ 
ting 

In this sub-section we focus toward TVA in its true sense. By bringing in an additional 
neighbor into consideration whenever v COM becomes zero, we propose an algorithmic way 
(Algorithm [1]) to implement TVA in a real system. Clearly, this provides a way to avoid 
ill-posedness associated with v COM being zero because non-zero speeds of individual agents 
ensure that considering an extra neighbor will make an otherwise zero v COM non-zero. More¬ 
over, this approach towards flocking can easily be restricted to a planar setting, and it has 
been demonstrated in what follows. 


5 Implementation Results 

5.1 Experimental Setup 

Our experimental test-bed is comprised of Pioneer 3 DX wheeled robots from Adept Mo- 
bileRobots |26j . These compact differential-drive mobile robots are equipped with reversible 
DC motors, high-resolution motion encoders and 19cm wheels, and the onboard computation 
is done via 32-bit Renesas SH2-7144 RISC microprocessor, including the P3-SH microcon¬ 
troller with ARCOS. The sensors on the robot include eight forward-facing ultrasonic (sonar) 
sensors. ARIA [27] provides an interface for controlling and receiving data from the robot, 
and communication with the robot for sending control commands (forward velocity and turn¬ 
ing rate ) is done via 802.11-b/g/n networking. The width of the robot is 380 mm and it has 
a swing radius of 260 mm. 

Algorithm implementation (i.e, feedback law computation) has been done in C++ using 
ROS [28], along with ROS-ARIA, as the interfacing robotics middleware. The experiments 
have been carried out in a laboratory environment equipped with a sub-millimeter accurate 
Vicon motion capture system ra- We use a Dell workstation to run ROS, and this computer 
is connected to the Vicon server via a dedicated Ethernet connection. 

The Vicon system captures the motion of the robots and sends out the position and 
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Algorithm 1: Topological Velocity Alignment - 3D 
Data: Initial Time - t imtm i ; Final Time - tfi na i ; Sampling Interval - A; Number of 
Agents - n\ Initial Position and Orientation - {<7i}” =1 ; Neighborhood Size - K 


begin 

Initialize. t curren t * tinitial > 

for i = 1 to n do 
| Initialize: State - A', <— c 


while t curren t $5 tfinai do 

for i = 1 to n do 

Define: M, - the set of A-nearest neighbors ; 

Compute: Neighborhood center of mass velocity - v COM ; 

if V COM = 0 then 

Define: J\T t - the set of K + 1-nearest neighbors ; 
Compute: Neighborhood center of mass velocity - v COM ; 

Compute: Steering control - Ui , Vi ; 

Implement: Steering Control - {ui,Vi}™ =1 ; 

Update: State - {Aj}” =1 ; 

Update. Time - t curren t < t curren t ~ I A . 


heading data to the computer running ROS. The control law program listens to this data, 
and transmits the individual turning rates (as individual speeds are assumed to be constant). 
Both of these operations are carried out at a frequency of 25 Hz. As the velocity vector (rf, 
with k denoting the time index) is aligned along the robot heading, x( : and yf can be directly 
computed from the heading data. Then, the curvature variable u k is evaluated from the 
corresponding control laws (JH1 CE2]), and the turning rate u> k {i = 1,2 ,,n) (in degrees/sec) 
is computed as: 


uj, 



(32) 


Next we will present the implementation results of the two control laws in our robotic 
test-bed. In this paper, we are presenting results for which the speeds of all the individual 
agents are same, i.e. ^ = uj,\/i,j. Though it should be noted that both control laws can be 
implemented with different speeds. 


5.2 Implementation of MMC 

Here we will describe the implementation results of MMC, difficulties of pure MMC law (|5il 
and usefulness of using the dissipative control term (j8lh First, the original MMC law (J5]) is 
used to control the robots. The observed trajectories are then compared with theoretically 
predicted ones. To generate the ideal trajectories in discrete time, we integrate the reduced 
system dynamics ([6]) . Considering the conserved quantity in the system (J7|), we used the 
method described in (30j for integration instead of general ODE solver, which otherwise would 
not be able to keep the quantity E(p, 7) constant to its initial value. From the integrated 
values of p k and 7 fc , we reconstruct the trajectories (i.e. r; fc , Xi fc , yi fc for i = 1,2) |l3j . where 
k denotes the time index. At each time instance t k , the error (ef) is then calculated as: 
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(a) Trajectory (Experiment) 


(b) Trajectory (Simulation) 



(c) Phase plot 


(d) Curvature 
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Figure 4: Comparison between experimental and simulation results under pure MMC law in 
action (with fx = 0.001 mm” 1 and Vi = v 2 = 200 mm/sec) 


pk _ _ yk I 

1 1 i,expt 1 i,ideal I * 

The plots of a sample run using pure MMC law (0 are given in FigHJ The parameters 
for this run are, /! = 0.001 mm” 1 and V\ = V 2 = 200 mm/sec. It can be seen from the figures 
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(a) Trajectory (Experiment) 



(c) Phase plot 
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(b) Trajectory (Simulation) 



(d) Curvature 



Time (in sec) 

(f) Percentage error in conserved quantity 



Figure 5: Performance of the system significantly improved upon addition of the dissipative 
control term (with /i = 0.001 mm -1 , U\ = u 2 = 200 nnn/sec and kd = 1 x 10~ 15 mm~ 6 sec 3 ) 

that the experimental trajectories diverge from the ideal ones and thus the error also keeps 
increasing over time in Fig 03 This also affects the phase plot in Fig 03 as we can clearly see 
the spiraling out type behavior as opposed to the ideal periodic behavior. Possible reasons 
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behind this have already been discussed in section 13.11 

To overcome the problem with original MMC law, a stabilization term has been intro¬ 
duced in the control law (151). We applied the modified control law with kd = 1 x 10~ 15 
mm~ 6 sec 3 and kept all other parameters fixed. The resulting performance is quite satisfac¬ 
tory as shown in Fig. [5a] (refer |3T] for implementation video). The error (Fig l5e|) is also 
bounded (~ 250 mm) within the size of the robots (~ 400 mm). Also, it can be seen from 
Fig [5f| that the cumulative average of the error in conserved quantity converges to almost 
0%, indicating superior performance. 

5.3 Implementation of TVA 

We implemented the TVA control law (115]) in a 2 dimensional setting (i.e. Vi(t) is ignored). 
As the implementation is in discrete time, we followed Algorithm [Tj in our implementation 
in order to avoid the singular case of |v COM | = 0. To demonstrate the performance of TVA 
control law, we designed three different experiments (refer [3T] for implementation videos). 
In these experiments, the sonar sensors on the robots were activated to sense any obstacle 
in the direction of motion of the robots and if any robot can sense such an obstacle, it will 
simply apply a maximum turning rate ( ui sat ) to avoid collision. The sonars are programmed 
to detect an obstacle only in close proximity (~ 300 mm) of the robots. In all our experiments 
co sat is taken to be 50 rad/sec, forward speeds of all of the robots are kept constant at 60 
mm/sec and the value of the parameter // is chosen to be 1 Hz. 

5.3.1 Experiment 1 

A system with eight agents is considered and we apply same TVA law to all of them. 
The neighborhood size is taken to be three (i.e. I\ =3). The robots are initially placed 
in arbitrary positions and directions. The footprints of the robots and the corresponding 
contrast function, 0(f) = is plotted against time in Fig [6al [60 The initial and 

final directions of the robots are shown using arrows and the final positions of the robots 
are denoted using dots. It can be seen from Fig [6b] that the contrast function decays to zero 
very quickly which indicates perfect velocity alignment within the swarm. 

5.3.2 Experiment 2 

Next we decreased the neighborhood size and made it one (K = 1), so that every robot 
only ‘communicates’ with its closest neighbor. We chose the initial positions in such a way 
that they may form sub-clusters instead of moving as a single swarm. This behavior is 
called ‘ splitting 1 in a swarm. From Fig [63 we can clearly see that the swarm of eight robots 
gradually split from each other and form three different clusters. It is to be noted that even if 
all the agents are not going in the same direction, the contrast function still converges to zero 
(Fig. l6dl) . This happens because each of the robots are aligned with their nearest neighbors 
and hence each of the individual contrast functions (0*(f)) are zero. This experiment may 
explain the splitting phenomenon observable in nature. 

5.3.3 Experiment 3 

Lastly, we combined the above two experiments, and conducted an experiment using six 
robots in a swarm and another robot as a predator. A separate computer was used for 
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(a) Trajectories: 8 agents, flocking 


(c) Trajectories: 8 agents, splitting 





Time (in sec) 


Time (in sec) 

(d) Contrast function: 8 agents, splitting 




(e) Trajectories: 6 agents, perturbation 


(f) Contrast function: 8 agents, perturbation 


Figure 6: Trajectory and contrast functions of TVA for (i) Experiment 1 (a),(b) with 8 
agents, demonstrates flocking behavior (K = 3); (ii) Experiment 2 (c),(d) with 8 agents, 
describes the splitting behavior due to low neighborhood size (K = 1), and (iii) Experiment 
3 (e),(f) with 6 agents, shows perturbation can cause a swarm to split, the trajectory of 
the perturbing agent is not shown. (// = 1 Hz and zzj = 60 mm/sec is kept fixed for all 
experiments) 
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manual control of the ‘predator’ robot. 

At the beginning, neighborhood size is kept at K = 3, such that the ‘communication’ 
graph among the robots stays connected and they move as an entire swarm in a common 
direction. When the swarm comes close to the predator, the neighborhood size is decreased 
to one. As we are not using any onboard visual sensing and the sonar sensing is done only 
in very close region (~ 300 mm), the change in neighborhood size is made manually. From 
Fig |6fJ we can see that the change in neighborhood size takes place at around 20 seconds 
and we can also see a tiny jump in the contrast function at that time. The predator then 
slowly approaches to one of the agents in the swarm, which abiding to its collision avoidance 
rule, turns to avoid the predator. In Fig [6^ the trajectories of the agents are drawn in 
dashed lines before the occurrence of this event and in solid lines afterwards. The trajectory 
of the predator robot in not shown in the figure. After creating the initial perturbation, 
the predator is slowly moved through the swarm causing some subsequent disturbances. 
These perturbations create a noticeable impact in the swarm. As the attacked agent turns, 
its neighbor also tries to align itself with that agent and so does its neighbor. This goes 
on until the communication graph becomes disconnected and a split in the swarm is then 
observed [3CE] just like in Experiment 2. As we can see in Fig [1x3 the swarm is divided in two 
groups after the attack of the predator. The jumps in plot of the contrast function in Fig 
l6fl symbolize the perturbations caused by the external agent. The contrast function value 
eventually converges to zero after the members are aligned with their neighbors within each 
subgroup. 

6 Conclusion and Future Work 

This paper has introduced a control strategy (TVA) which, based on local information, at¬ 
tempts to align the individual velocities, and the theoretical analysis for a special case of 
planar two-agent system has been complemented by experiments in a laboratory environ¬ 
ment. Additionally, we chose MMC for implementation due to its dynamic coverage property 
with potential in many domains (resource harvesting, environmental monitoring, search and 
rescue missions etc.). Moreover, the exploitive nature of MMC (through coverage in a certain 
annular region) augments the exploratory nature of TVA (through making the agents move 
in a common direction). Future works will primarily focus on two areas, namely effect of a 
beacon (influencing both agents) into MMC, and understanding the effect of a covert leader 
(equipped with extra information) in a flock. Also, this study does not address how the 
system behavior will change in a noisy (or uncertain) environment. In future, the authors 
have plans to extend current results to study the effect of sensor noise (along the lines of 
[32]) and perform robustness analysis. 
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